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ABSTRACT 

The  necessary  and  sufficient  conditions  in  terms  of  antenna  array  geometry  for  unique 
estimation  of  azimuth-elevation-carriers  (AEC)  of  narrow-band  signals  are  introduced. 
These  conditions  make  it  possible  to  solve  the  problem  of  joint  AEC  estimation  based 
on  spatial  samples  taken  by  the  volume  array  without  using  the  large  number  of 
temporal  samples. 

INTRODUCTION 

The  problem  of  joint  azimuth-elevation-carriers  (AEC)  estimation  of  narrow-band 
signals  arises  in  communication,  radar,  radio  astronomy,  sonar  and  many  other 
applications.  An  important  issue  related  to  this  problem  is  the  unique  estimation  of 
signal  AECs  via  array  of  given  geometry.  This  issue  called  also  as  the  identifiability 
problem  (see  [l]-[3]  and  references  herein).  The  standard  formulation  of  this  problem 
involves  determination  of  the  number  of  harmonics  that  can  be  resolved  for  a given  total 
sample  size.  The  solutions  to  be  proposed  are  obtained  under  conditions  of  equispaccd 
samples,  existence  of  all  samples  along  all  the  array  dimensions  and  difference  of  all 
harmonics  to  be  resolved.  These  conditions  significantly  restrict  the  practical 
implementation  of  the  theory.  In  this  paper  we  present  the  necessary  and  sufficient 
conditions  for  unique  AEC  estimation  by  volume  array  those  are  given  in  terms  of  the 
number  of  the  antenna  array  sensors  and  array  geometry. 

DATA  MODEL  AND  PROBLEM  FORMULATION 

We  assume  that  there  are  M point  sources  that  are  emitting  the  unknown  narrow-band 
complex  deterministic  signals  sm(t),  in  the  direction  of  the  measurement 

system.  For  every  signal  the  parameters  of  interest  are  the  azimuth  p„, , the  elevation 

s,„  and  the  carrier  /,,.  To  describe  the  signal  parameters  let  us  introduce  a vector 

Mm  = fmem  If  max  G R?' - Mm  e U » = {m  - P e [O.  71 J £ G [O.  7t  J / G , /„„„.  ]}  , where 

e,„=[cosfimco.sE„h  .v/7?P„,c'o.v£„„  sin(:n,]!  is  the  unit  vector  in  the  Cartesian  coordinate 
system  and  (-)7  denotes  transpose.  It  is  assumed  that  the  measurement  system  consists  of 
N omnidirectional  point  sensors.  Then  the  model  of  antenna  array  output  in  noiseless 
case  can  be  represented  as  follows: 

z(/)  = Ts(/)eCAxl,  / = 1...L.  (1) 
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where  A=[a(\i.i)...  «(li«/)|^(  ' is  a matrix  of  steering  vectors,  a(p,„)=[l,  exp {]2n 

d2\im},  •••,  exp{j27i  dj\xm}]'  is  a steering  vector,  dn  = dninc / fmax  is  a vector  of  the  nth 

sensor  position  ( d]=0 ),  dn  is  a distance  to  the  nth  sensor,  in  is  a unit  vector  toward  the 

nth  sensor  position,  D=[di„.  rfy] e R3xN  is  a matrix  that  specify  the  antenna  array 
geometry,  s(t)eCMxI  is  a vector  of  the  signal  waveforms,  L denotes  the  number  of  data 
snapshots  available.  We  assume,  that  the  sampled  covariance  matrix  of  signal 

waveforms  P = (L)~xlL^xs{t)sH  {t)<=.  CM'M  has  full  rank,  that  is  rank(p)=  M , and  the 
number  of  signals  M is  known. 

The  problem  to  be  solved  is  formulated  as  follows:  to  determine  the  basic  requirements 
for  the  number  N and  positions  D of  array  sensors,  that  guarantee  the  unique  estimation 
of  any  distinct  set  /q  * ...  * fiM  e U of  signal  parameters. 

IDENTIFIABILITY  OF  AEC  ESTIMATES 

The  necessary  and  sufficient  conditions  for  identiflability  of  U are  [1]  NSCN : the  set 
of  steering  vectors  «(p)  is  known;  NSC2:  for  any  fix*...*nM€.U  the  matrix  A has 

full  rank  rank[a(/ix)...a(fiM  )]=  min(N, M)  = M < N ; NSC3:  the  number  of  sensors  N, 
the  rank  R of  matrix  P and  the  number  Q of  parameters  per  signal  satisfy  to  the 
following  inequality  N > R + Q = M + 3.  The  NSC]  and  NSC3  hold  if  the  antenna  array  is 
calibrated  and  its  geometry  D is  known,  and  the  maximum  number  of  impinging  signals 
M < N - 3 . Therefore,  the  problem  to  be  solved  can  be  reduced  to  the  determination  of 
the  requirements  for  the  sensor  positions  D that  guarantee  the  implementation  of  the 
NSC2. 

Theorem : For  any  /q  * ...  * fiM  e Uyi  the  matrix  A has  full  rank  if:  i)  the  array  has  a 
subarray  with  K <N  sensors,  such  that  the  number  N p of  the  parallel  planes  that  can 
be  passed  through  the  points  dv...,dK  satisfy  the  conditions  Np  = K -2>  M , and  ii) 
the  maximum  distance  dmax  between  any  nearest  sensors  of  subarray  is  less  then  half  of 
minimum  wavelength  dmax  < Xmm  / 2 = c/{2fmax). 

Proof.  The  proof  is  based  on  the  fact  that  the  rank  of  matrix  A is  defined  as  follows 
rank(A)  = min(N  - Kr,M  - Kc),  where  Kr  is  the  number  of  rows  and  Kc  is  the 

number  of  columns  that  are  coincident  to  any  other  row  or  column  respectively. 
Condition  ii)  means,  that  the  maximum  sampling  frequency  in  spatial  domain  is  higher 
then  the  Nyquist  rate,  therefore  the  vectors  p giving  Kc  > 0 columns  that  are 

coincident  to  the  column  fl(/q ) are  the  roots  of  the  following  set  of  equations 

iTn]i  = cn,  n = 2...K , (2) 

where  cn  = /J/q  . The  vectors  /r"1  that  satisfy  to  the  i?th  equation  of  (2)  represent  the 
set  of  points  lying  on  a plane  z(in,cn)eR2  that  is  orthogonal  to  the  vector  in  and 
located  on  the  distances  \cn  | < 1 from  the  origin.  Therefore  the  set  U’  of  roots  of  (2) 
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are  the  intersection  between  U and  the  planes  s(in,c„): 
U'Ji 2-E)=U„ni2E(/,„c„). 

Observe,  that  the  set  U is  a space  between  two  hemispheres  of  radiuses  fmin/fm„x 
and  1.  Therefore,  if  the  vectors  dx...dK  do  not  lie  on  the  same  plane,  i.e. 
ranked x ...d K ])  = 3 , then  the  set  U'  is  a point  and  only  point  U'  = {//, } . It  means,  that 
if  the  subarray  is  a volume,  then  for  any  px  *...■&  ft Xl  e(/p  the  inequality 
«(//, ) ^ ...  * ) holds  true  and  Kc  = 0 . Observe  also,  that  if  < K sensors  lie  on 

a plane,  then  the  set  U'^  [i2...iK^ ) is  a section  and  consists  of  infinite  number  of 
elements.  Hence,  for  any  distinct  ^...^//A/  e U'^  [i2...iK  ),  A,  equations  of  (2) 

hold  and  K}  rows  of  the  matrix  [«(//, )...«(// A/ )]  are  equal,  i.e.  Kr  = K ,.  The  K r -1 
rows  can  be  rejected  without  loss  of  information.  Therefore,  to  ensure  the  full  rank  of 
the  matrix  A the  number  of  sensors  lying  in  nonparallel  planes  has  to  satisfy  the 
condition  N > M . Since  the  plane  can  be  passed  through  any  three  sensors,  the 

minimum  size  of  the  subarray  is  K = N +2  . 

CONCLUSION 

An  important  corollary  of  the  theorem  is  that  the  joint  AEC  estimation  of  narrow-band 
signals  can  be  performed  by  means  of  spatial  samples  taken  by  volume  antenna  array. 
The  temporal  averaging  is  needed  to  ensure  a full  rank  of  the  covariance  matrix  of 
signal  waveforms  as  well  as  to  increase  the  signal-to-noise  ratio.  In  practice  analogue 
tools  can  successfully  carry  out  this  procedure  that  is  very  important  for  the  MW 
measurements  systems.  Another  conclusion  is  that  the  arrays  with  irregular  geometry 
have  advantages  in  comparison  with  the  periodic  arrays  from  the  identifiability  point  of 
view. 

The  possible  applications  of  the  considered  theory  could  include  the  wide-range 
spectrum  monitoring  direction  finders,  Doppler  radar  systems,  wireless  communication 
systems  and  other. 
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